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ABSTRACT 


Recently,  Goldberg  proposed  a  new  approach,  to  the  maximum  network 
^  flow  problem.  The  approach  yields  a  very  simple  algorithm  running  in  0(n3) 
‘'  time  on  n-vertex  networks.  Incorporation  of  the  dynamic  tree  data  structure  of 
Sleator  and  Tarjan  yields  a  more  complicated  algorithm  with  a  running  time  of 
0(rtm  log  ( n2/m ))  on  m-e dge  networks/ Ahuja  and  Orlin  developed  a  variant  of 
Goldberg’s  algorithm,  that  uses  scaling  and  runs  in  0(nm  +  n2  log U)  time  on  net¬ 
works  with  integer  edge  capacities  bounded  by  U.  ty  this  paper  we-first  obtairv'a 
modification  of  the  Ahuja-Orlin  algorithm  with  a  running  time  of 

0(nm  +  n2-~^~).  We  then-show  that  the  use  of  dynamic  trees  in  this  algo- 
loglogt/ 

rithm  further  reduces  the  time  bound  to  0(nm  log(— -n  +  2)).  This  result 

m  loglogt/ 

demonstrates  that  the  combined  use  of  scaling  and  dynamic  trees  results  in  speed 
not  obtained  by  using  either  technique  alone. 
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Improved  Time  Bounds  for  the  Maximum  Flow  Problem 


Ravindra  K.  Abuja  1,2 
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1.  Introduction 

We  consider  algorithms  for  the  classical  maximum  network  flow  problem  [4,5,11,13,18]. 
We  formulate  the  problem  as  follows.  Let  G  =  (V,E)  be  a  directed  graph  with  vertex  set  V  and 
edge  set  £.  The  graph  G  is  a  flow  network  if  it  has  two  distinct  distinguished  vertices,  a  sink  s  and 
a  source  t,  and  a  non-negative  real-valued  capacity  u(y,w )  on  each  edge  (v,w)  e  E.  We  assume 
that  G  is  symmetric,  i.e.  (v,w)  e  E  iff  (w,v)  e  £.  We  denote  by  n,m,  and  U  the  number  of  ver¬ 
tices,  the  number  of  edges,  and  the  maximum  edge  capacity,  respectively.  For  ease  in  stating 
time  bounds,  we  assume  m  t  n- 1  and  U  I>  4.  Bounds  containing  U  are  subject  to  the  assumption 
that  all  edge  capacities  are  integral.  All  logarithms  in  the  paper  are  base  two  unless  an  explicit 
base  is  given. 

A  flow  f  on  a  network  G  is  a  real-valued  function  /  on  the  edges  satisfying  the  following 
constraints: 


/(v,w)  £  u(v,w)  for  all  (v.w)  e£  (capacity  constraint),  (1) 

f(v,w)  =  - f(w,v )  for  all  (v.w)  e  £  (antisymmetry  constraint),  (2) 

£  /(v.w)  =  0  for  all  w  e  V-{s,f )  (conservation  constraint).  (3) 

v:(v,w)cE 
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The  value  I /I  of  a  flow  /is  the  net  flow  into  the  sink: 


1/1=  I 

»:(»,<)  e  E 


A  maximum  flow  is  a  flow  of  maximum  value.  The  maximum  flow  problem  is  that  of  finding  a 
maximum  flow  in  a  given  network. 

Remark.  We  assume  that  all  edge  capacities  are  finite.  If  some  edge  capacities  are  infinite  but  no 
path  of  infinite-capacity  edges  from  s  to  t  exists,  then  each  infinite  capacity  can  be  replaced  by  the 
sum  of  the  finite  capacities,  without  affecting  the  problem.  □ 

The  maximum  flow  problem  has  a  long  and  rich  history,  and  a  series  of  faster  and  faster 
algorithms  for  the  problem  have  been  developed.  (See  [9]  for  a  brief  survey.)  The  previously 
fastest  known  algorithms  are  those  of  Goldberg  and  Taijan  [7,9],  with  a  running  time  of 
0{nm  log  ( n2lm )),  and  that  of  Ahuja  and  Orlin  [1],  with  a  running  time  of  0(nm  +  n2  lo gU). 
Both  of  these  algorithms  are  refinements  of  a  generic  method  proposed  by  Goldberg  [7],  which 
we  shall  call  the  preflow  algorithm.  For  networks  with  m  =  Q(n2),  the  Goldberg-Tatjan  bound  is 
0(n3),  which  matches  the  bound  of  several  earlier  algorithms  [10,12,14,19].  For  networks  with 
m  =  0(n2-c)  for  some  constant  e  >  0,  the  Goldberg-Taijan  bound  is  0(nm  logn),  which  matches 
the  bound  of  the  earlier  Sleator-Taijan  algorithm  [15,16].  Under  the  similarity  assumption  [6], 
namely  U  =  0(n  *)  for  some  constant  k,  the  Ahuja-Oriin  bound  beats  the  Goldberg-Tatjan  bound 
unless  m  -  0(n)  or  m  =  Q(n2). 

The  Goldberg-Taijan  and  Ahuja-Oriin  algorithms  obtain  their  speed  from  two  different 
techniques.  The  former  uses  a  sophisticated  data  structure,  the  dynamic  tree  structure  of  Sleator 
and  Taijan  [16,17,18],  whereas  the  latter  uses  scaling.  Our  main  purpose  in  this  paper  is  to 
demonstrate  that  the  use  of  both  techniques  results  in  efficiency  not  obtained  by  using  either  tech¬ 
nique  alone.  We  first  modify  the  Ahuja-Oriin  algorithm  to  obtain  a  small  improvement  in  running 

time,  to  0(nm  +  n2  )  Then  we  show  that  the  use  of  dynamic  trees  in  the  modified  algo- 

loglogt/ 

rithm  results  in  a  running  time  of  0(nm  Iog( — *1  +  2)).  Under  the  similarity  assumption, 

m  loglogU 

this  bound  is  better  than  all  previously  known  bounds  for  the  maximum  flow  problem.  For  net¬ 
works  with  m  =  O(n)  and  U  =  0(nk)  for  some  constant  k,  the  bound  is  0(nm  log  logn),  which 
beats  both  the  Goldberg-Taijan  bound  and  the  original  Ahuja-Oriin  bound  by  a  factor  of 
logn/loglogn.  Moreover,  the  bound  rapidly  approaches  0(nm )  as  the  graph  density  m/n 
increases. 


Our  paper  consists  of  four  sections  in  addition  to  this  introduction.  In  Section  2  we  review 
the  preflow  algorithm.  In  Section  3  we  review  the  Ahuja-Orlin  algorithm  and  describe  and 
analyze  our  improvement  of  it.  In  Section  4  we  combine  the  use  of  dynamic  trees  with  the 
method  of  Section  3  and  analyze  the  resulting  algorithm.  Section  5  contains  some  final  remarks. 

2.  The  Preflow  Algorithm 

In  contrast  to  the  classical  augmenting  path  method  of  Ford  and  Fulkerson  [5],  which  moves 
flow  along  an  entire  path  from  s  to  t  at  once,  the  prefiow  method  moves  flow  along  a  single  edge 
at  a  time.  The  key  concept  underlying  the  algorithm  is  that  of  a  preflow,  introduced  by  Karzanov 
[10].  A  preflow  f  is  a  real-valued  function  on  the  edges  satisfying  constants  (1),  (2),  and  a  relaxa¬ 
tion  of  (3).  For  any  vertex  w,  let  the  flow  excess  of  w  bee  (w)  =  £  /(v,w). 

v:(v,H>)e  E 

The  required  constraint  is  the  following: 

e(w)  >  0  for  all  w  e  V-[ s }  (nonnegativity  constraint).  (4) 

We  call  a  vertex  v  active  if  v*  r  and  e(y)  >  0.  Observe  that  the  nonnegativity  constraint  implies 
that  e(s)  <  0. 

The  residual  capacity  of  an  edge  (y,w)  with  respect  to  a  preflow  /  is 
ufly.w)  =  u(y,w)-f(v,w).  An  edge  is  saturated  if  ufly.w)  =  0  and  unsaturated  otherwise.  (The 
capacity  constraint  implies  that  any  unsaturated  edge  (v.w)  has  ufly.w)  >  0.) 

The  preflow  algorithm  maintains  a  preflow  and  moves  flow  from  active  vertices  through 
unsaturated  edges  toward  the  sink,  along  paths  estimated  to  contain  as  few  edges  as  possible. 
Excess  flow  that  cannot  be  moved  to  the  sink  is  returned  to  the  source,  also  along  estimated  shor¬ 
test  paths.  Eventually  the  preflow  becomes  a  flow,  which  is  a  maximum  flow. 

As  an  estimate  of  path  lengths,  the  algorithm  uses  a  valid  labeling,  which  is  a  function  d 
from  the  vertices  to  the  nonnegative  integers  such  that  d(s)  =  n,  d(t)  =  0,  and  d(y)  S  d(w)  +  1  for 
every  unsaturated  edge  (v,w).  A  proof  by  induction  shows  that,  for  any  valid  labeling  d. 
d(y)  $  min  {dflv.s)  +  n.dflv.t)},  where  dflv.w )  is  the  minimum  number  of  edges  on  a  path  from 
v  to  w  consisting  of  edges  unsaturated  with  respect  to  the  flow/.  We  call  an  edge  (v,h>)  eligible  if 
(v,h*>  is  unsaturated  and  d(v)  =  d(w)  +  1. 


Hie  algorithm  begins  with  an  initial  preflow  f  and  valid  labeling  d  defined  as  follows: 

u(v,w)  if  v=s, 
f(v,w)  =•  -u(w,v)  if  w  =  s, 

0  if  v  *  s  and  w  *  s. 


d(y)  =  min  [df  (v,s)  +  n,dj(v,t)) . 

The  algorithm  consists  of  repeating  the  following  two  steps,  in  any  order,  until  no  vertex  is  ac¬ 
tive: 

push  (v.w). 

Applicability:  Vertex  v  is  active  and  edge  (v,w)  is  eligible. 

Action:  Increase  f(v,w )  by  min  {e(v),uy{v,w)}.  The  push  is  saturating  if  (v,w)  is 
saturated  after  the  push  and  nonsaturating  otherwise. 

relabel  (v). 

Applicability:  Vertex  v  is  active  and  no  edge  (v,w)  is  eligible. 

Action:  Replace  d(v)  by  min  [d(w)  +1  I  (v,w)  is  unsaturated}. 

When  the  algorithm  terminates,  / is  a  maximum  flow.  Goldberg  and  Taijan  derived  the  fol¬ 
lowing  bounds  on  the  number  of  steps  required  by  the  algorithm: 

Lemma  2.1  [9].  Relabeling  a  vertex  v  strictly  increases  d(y).  No  vertex  label  exceeds  2n-l.  The 
total  number  of  relabelings  is  0(n2). 

Lemma  2.2  [9].  There  are  at  most  0(nm)  saturating  pushes  and  at  most  0(n2m)  nonsaturating 
pushes. 

Efficient  implementations  of  the  above  algorithm  require  a  mechanism  for  selecting  pushing 
and  relabeling  steps  to  perform.  Goldberg  and  Taijan  proposed  the  following  method.  For  each 
vertex,  construct  a  (fixed)  list  A(v)  of  the  edges  out  of  v.  Designate  one  of  these  edges,  initially 
the  first  on  the  list,  as  the  current  edge  out  of  v.  To  execute  the  algorithm,  repeat  the  following 
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step  until  there  no  active  veitices: 

push/relabel  (v). 

Applicability:  Vertex  v  is  active. 

Action:  If  the  current  edge  (v,w)  of  v  is  eligible,  perform  push(y,w).  Otherwise,  if  (v,w) 
is  not  the  last  edge  on  A(v),  make  the  next  edge  after  (v,w)  the  current  one.  Other¬ 
wise,  perform  relabel  (v)  and  make  the  first  edge  on  A(v)  the  current  one. 

With  this  implementation,  the  algorithm  runs  in  0(nm )  time  plus  0(1)  time  per  nonsaturat¬ 
ing  push.  This  gives  an  0(n2m)  time  bound  for  any  order  of  selecting  vertices  for  push/relabel 
steps.  Making  the  algorithm  faster  requires  reducing  the  time  spent  on  nonsaturating  pushes.  The 
number  of  such  pushes  can  be  reduced  by  selecting  vertices  for  push/relabel  steps  carefully. 
Goldberg  and  Taijan  showed  that  FIFO  selection  (first  active,  first  selected)  reduces  the  number 
of  nonsaturating  pushes  to  0(n3).  Cheriyan  and  Maheshwari  [2]  showed  that  highest  label  selec¬ 
tion  (always  pushing  flow  from  a  vertex  with  highest  label)  reduces  the  number  of  nonsaturating 
pushes  to  0(n2mv2).  Ahuja  and  Orlin  proposed  a  third  selection  rule,  which  we  discuss  in  the 
next  sectioa 

3.  The  Scaling  Preflow  Algorithm 

The  intuitive  idea  behind  the  Ahuja-Oriin  algorithm,  henceforth  called  the  scaling  preflow 
algorithm,  is  to  move  large  amounts  of  flow  when  possible.  The  same  idea  is  behind  the  max¬ 
imum  capacity  augmenting  path  method  of  Edmonds  and  Karp  [3]  and  the  capacity  scaling  algo¬ 
rithm  of  Gabow  [6].  One  way  to  apply  this  idea  to  the  preflow  algorithm  is  to  always  push  flow 
from  a  vertex  of  large  excess  to  a  vertex  of  small  excess,  or  to  the  sink.  The  effect  of  this  is  to 
reduce  the  maximum  excess  at  a  rapid  rate. 

Making  this  method  precise  requires  specifying  when  an  excess  is  large  and  when  small. 
For  this  purpose  the  scaling  preflow  algorithm  uses  an  excess  bound  A  and  an  integer  scaling  fac¬ 
tor  k  t  2.  A  vertex  v  is  said  to  have  large  excess  if  its  excess  exceeds  A  Ik  and  small  excess  other¬ 
wise.  As  the  algorithm  proceeds,  k  remains  fixed,  but  A  periodically  decreases.  Initially,  A  is  the 
smallest  power  of  k  such  that  At  U.  The  algorithm  maintains  the  invariant  that  e(v)  <,  A  for  every 
active  vertex  v.  This  requires  changing  the  pushing  step  to  the  following. 


push  (v.w). 

Applicability:  Vertex  v  is  active  and  edge  (v,w)  is  eligible. 

Action:  If  w  *  r,  increase /(v,w)  by  min  (e(v),u/  (v,w),  A-e(w)}.  Otherwise  (w  =  t). 
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incnease/(v,H')by  min  { e(v),Uf(v,w)}. 

The  algorithm  consists  of  a  number  of  scaling  phases,  during  each  of  which  A  remains  con¬ 
stant.  A  phase  consists  of  repeating  push/relabel  steps,  using  the  following  selection  rule,  until 
no  active  vertex  has  large  excess,  and  then  replacing  A  by  A/k. 

Large  excess,  smallest  label  selection:  Apply  a  push/relabel  step  to  a  vertex  v  of  large  excess; 
among  such  vertices,  choose  one  of  smallest  label. 

If  the  edge  capacities  are  integers,  the  algorithm  terminates  after  at  most  [log*  U  +  lj 
phases:  after  [log*  U  +  lj  phases,  A  <  1,  which  implies  that /is  a  flow,  since  the  algorithm  main¬ 
tains  integrality  of  excesses.  Ahuja  and  Orlin  derived  a  bound  of  0(kn2  log*  U)  on  the  total 
number  of  nonsaturating  pushes.  We  repeat  the  analysis  here,  since  it  provides  motivation  for  our 
modification  of  the  algorithm. 

Lemma  3.1  [1].  The  total  number  of  nonsaturating  pushes  in  the  scaling  preflow  algorithm  is 
0(kn2  (log*  U  +  1)). 

Proof.  Consider  the  function  O  =  £  e(y)  d(v)l A.  We  call  <t>  the  potential  of  the  current 

v  active 

preflow  /  and  labeling  d.  Since  0  <  e(v)/A  <,  1  and  0  £  d(v)  <,  2 n  for  every  active  vertex  v, 
0  <,  <I>  <*  2n  throughout  the  algorithm.  Every  pushing  step  decreases  <I>.  A  nonsaturating  push¬ 
ing  step  decreases  <I>  by  at  least  1  Ik.,  since  the  push  is  from  a  vertex  v  with  excess  more  than  A/k 
to  a  vertex  w  with  d(w)  =  d(v)  -1 ,  and  e(w)£  A/k  or  w  =  t.  The  value  of  d>  can  increase  only  dur¬ 
ing  a  relabeling  or  when  A  changes.  A  relabeling  of  a  vertex  v  increases  d>  by  at  most  the  amount 
d(y)  increases.  Thus  the  total  increase  in  A  due  to  relabelings,  over  the  entire  algorithm,  is  at 
most  2« 2.  When  A  changes,  increases  by  a  factor  of  k,  to  at  most  2 n2.  This  happens  at  most 
[log*  U  +  lj  times.  Thus  the  total  increase  in  <D  over  the  entire  algorithm  is  at  most  2 n2 
l log*  U  +  2| .  The  total  number  of  nonsaturating  pushes  is  at  most  k  times  the  sum  of  the  initial 
value  of  O  and  the  total  increase  in  <I>.  This  is  at  most  2kn 2  [_  log*  U  +  3j  .  □ 

Choosing  it  to  be  constant  independent  of  n  gives  a  total  time  bound  of  0(nm  +  n2  log U) 
for  the  scaling  preflow  algorithm,  given  an  efficient  implementation  of  the  vertex  selection  rule. 
One  way  to  implement  the  rule  is  to  maintain  an  array  of  sets  indexed  by  vertex  label,  each  set 
containing  all  large  excess  vertices  with  the  corresponding  label,  and  to  maintain  a  pointer  to  the 
nonempty  set  of  smallest  index.  The  total  time  needed  to  maintain  this  structure  is 
0(nm  +  n2  log U). 
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Remark.  The  bound  for  the  scaling  preflow  algorithm  can  be  improved  slightly  if  it  is  measured 
in  terms  of  a  different  parameter.  Let  U*  =4  +  £  u  (s,v)/n.  Then  the  bound  on  nonsa- 

v:(j,v)e  E 

turating  pushes  can  be  reduced  to  O  (n2  log U*)  and  the  overall  time  bound  to 
O  ( nm  +  rt2  logi/*).  This  improvement  is  analogous  to  the  bound  Edmonds  and  Karp  derived  for 
their  capacity-scaling  transportation  algorithm  [3].  The  argument  is  as  follows.  The  algorithm 
maintains  the  invariant  that  the  total  excess  on  active  vertices  is  at  most  nU* .  Let  phase  i  be  the 
first  phase  such  that  A  £  U* .  Then  the  total  increase  in  <I>  due  to  phase  changes  up  to  and  includ- 

i 

ing  the  change  from  phase  i  -  1  to  phase  i  is  at  most  n  £  2n/2‘--'  =  0(n2).  The  total  increase  in 

h o 

O  due  to  later  phase  changes  is  O  ( n 2  log I/*).  □ 

Having  described  the  Ahuja-Orlin  algorithm,  we  consider  the  question  of  whether  its  run¬ 
ning  time  can  be  improved  by  reducing  the  number  of  nonsaturating  pushes.  The  proof  of 
Lemma  3.1  bounds  the  number  of  nonsaturating  pushes  by  estimating  the  total  increase  in  the 
potential  O.  Observe  that  there  is  an  imbalance  in  this  estimate:  0(n2  log*  U)  of  the  increase  is 
due  to  phase  changes,  whereas  only  0(n2)  is  due  to  relabelings.  Our  plan  is  to  improve  this  esti¬ 
mate  by  decreasing  the  contribution  of  the  phase  changes,  at  the  cost  of  increasing  the  contribu¬ 
tion  of  the  relabelings.  Making  this  plan  work  requires  changing  the  algorithm. 

We  use  a  larger  scale  factor  and  a  slightly  more  elaborate  method  of  vertex  selection.  We 
push  flow  from  a  vertex  of  large  excess,  choosing  among  such  vertices  one  with  highest  label. 
This  vertex  selection  rule  does  not  guarantee  that  a  nonsaturating  push  moves  enough  flow,  and 
we  need  a  stack-based  implementation  to  ensure  this.  The  algorithm  maintains  a  stack  S  of  ver¬ 
tices,  with  the  property  that  if  w  is  just  on  top  of  v  on  the  stack,  then  (v.w)  is  the  current  edge  out 
of  v.  At  the  beginning  of  a  scaling  phase,  the  stack  is  initialized  to  be  empty.  The  phase  consists 
of  repeating  the  following  step  until  no  vertex  has  large  excess,  and  then  replacing  A  by  A  Ik. 

stack  push! relabel.  If  S  is  empty,  push  onto  S  any  large-excess  active  vertex  of  largest  label. 

Now  S  is  definitely  nonempty.  Let  v  be  the  top  vertex  on  S  and  (v.hO  the  current  edge  out  of 
v.  Apply  the  appropriate  one  of  the  following  three  cases: 

Case  1:  (v,w)  is  eligible.  Perform  push  (v,w)  (modified  as  described  at  the  beginning  of 
this  section).  If  the  push  is  nonsaturating  and  e(v)  >  0,  push  w  onto  S.  Otherwise, 
while  S  is  nonempty  and  its  top  vertex  has  small  excess,  pop  S. 

Case  2:  (v,w)  is  not  eligible  and  is  not  last  on  A  (v).  Replace  (v,w)  as  the  current  edge 
out  of  v  by  the  next  edge  on  A  (v). 
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Case  3:  (v,w)  is  not  eligible  and  is  last  on  A  (v).  Relabel  v  and  pop  it  from  S.  Replace 
(v.wO  as  the  current  edge  out  of  v  by  the  first  edge  on  A(v).  While  S  is  nonempty  and 
its  top  vertex  has  small  excess,  pop  S. 

Observe  that  the  stack  push/relabel  step  is  just  a  push/relabel  step  augmented  with  some 
manipulation  of  the  stacks. 

Before  carrying  out  a  detailed  analysis  of  this  algorithm,  we  make  several  observations. 
Every  push  made  by  the  algorithm  is  from  a  vertex  of  large  excess.  A  push  along  an  edge  (v,w) 
results  in  one  of  three  events:  saturation  of  (v,w),  conversion  of  v  into  an  inactive  vertex,  or  addi¬ 
tion  of  w  to  5  with  an  excess  of  A.  In  the  third  case  (and  only  in  the  third  case),  the  push  may 
move  zero  flow. 

Lemma  3.2.  The  total  number  of  nonsaturating  pushes  made  by  the  stack-based  scaling  preflow 
algorithm  is  0(nm  +  kn 2  +  n2  Gog*  !/  +  1)). 

Proof.  To  bound  the  number  of  nonsaturating  pushes,  we  use  an  argument  like  the  proof  of 
Lemma  3.1,  but  with  two  potentials  instead  of  one.  The  first  potential  is  that  of  Lemma  3.1, 
=  £  «(v)  d(v)/A.  By  the  analysis  in  the  proof  of  Lemma  3.1,  every  push  decreases  <D.  and 

v  active 

the  total  increase  in  <I>  over  all  phases  is  0(n2  log*!/).  Every  nonsaturating  push  either  increases 
or  decreases  the  size  of  S.  Consider  a  nonsaturating  push  that  increases  the  size  of  S,  say  by  push¬ 
ing  flow  from  v  to  w  and  stacking  w.  When  w  is  stacked,  it  has  an  excess  of  A.  It  cannot  be 
unstacked  until  it  is  relabeled  or  until  its  excess  decreases  to  at  most  A  tk.  In  the  former  case  we 
charge  the  push  from  v  to  w  to  the  relabeling  of  w.  in  the  latter,  to  the  ai  least  1-1  Ik  decrease  in  d> 
caused  by  pushes  from  w  after  w  is  stacked  and  before  it  is  unstacked.  Since  1  —1  /A  ^  1/2  for  all 
fe  2,  there  are  a  total  of  0(n2)  pushes  charged  to  relabelings  and  0(n2  log *!/)  pushes  charged  to 
decreases  in  <I>. 

We  now  count  the  nonsaturating  pushes  that  reduce  the  size  of  S.  The  number  of  such 
pushes  that  begin  with  more  than  one  vertex  S  is  0(n2  log*!/),  since  each  such  push  must  be  pre¬ 
ceded  by  a  nonsaturating  push  that  adds  a  vertex  to  S,  and  the  number  of  these  is  0(n2  log*!/)  by 
the  above  argument. 

The  remaining  pushes  that  must  be  counted  are  those  that  begin  with  one  vertex  on  S  and 
end  with  none  on  S.  We  call  such  pushes  emptying.  An  emptying  push  can  cause  a  rather  small 
decrease  in  <1>,  namely  1/k;  thus  using  <I>  to  bound  such  pushes  only  gives  a  bound  of 
0(kn2  log*!/)  on  their  number.  We  count  emptying  pushes  more  carefully  by  using  a  second 
potential,  <t>2.  The  definition  of  <J>2  involves  two  parameters,  an  integer  l  and  a  set  P.  The  value 
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of  /  is  equal  to  the  minimum  of  2n  and  the  smallest  label  of  a  vertex  added  to  S  while  S  was 
empty  during  the  current  phase.  Observe  that  /  =  2n  at  the  beginning  of  a  phase  and  l  is  nonin¬ 
creasing  during  a  phase.  The  set  P  consists  of  all  vertices  with  label  greater  than  l  and  all  vertices 
with  label  equal  to  /  from  which  an  emptying  push  has  been  made  during  the  current  phase. 
Observe  that  P  is  empty  at  the  beginning  of  a  phase  and  P  never  loses  a  vertex  during  a  phase. 

The  definition  of  <t>2  is 


e(v)  (d(v)  -l  +  1)/A.  (IfP  =  0,<I>2=O.) 


V  E  P  :e(v)  >  0 


*1 


Observe  that  0  £  <t>2  ^  2n2.  Any  emptying  push  either  adds  a  vertex  to  P  or  decreases  d>2 
by  at  least  1  Ik.  The  number  of  vertices  added  to  P  is  at  most  n  [log*  U  +  lj  over  all  phases,  and 
hence  so  is  the  number  of  emptying  pushes  that  do  not  decrease  <l>2  by  at  least  1  Ik. 

To  bound  the  number  of  emptying  pushes  that  decrease  <h2  by  at  least  1/A:,  we  bound  the 
total  increase  in  C>2.  Increases  in  <h2  are  due  to  relabelings  and  to  decreases  m  l.  (A  vertex  added 
to  P  because  of  an  emptying  push  has  zero  excess  and  hence  adds  nothing  to  <J>2.)  A  relabeling  of 
a  vertex  v  increases  <I>2  by  at  most  the  increase  in  d(v)  plus  one;  the  ‘‘plus  one”  accounts  for  the 
fact  that  the  relabeling  may  add  v  to  P.  Thus  relabelings  contribute  at  most  4 n2  to  the  growth  of 
<b2. 

There  are  at  most  In  decreases  in  l  per  phase.  A  decrease  in  /  by  one  adds  at  most  nlk  to 
<h2,  since  when  the  decrease  occurs  every  vertex  in  P  has  small  excess.  Thus  the  total  increase  in 
<I>2  due  to  decreases  in  l  is  at  most  In 2  [log*  U  +  lj  Ik  over  all  phases. 

The  total  number  of  emptying  pushes  that  decrease  <I>2  by  at  least  1/A:  is  at  most  k  times  the 
total  increase  in  <X>2,  since  <I>2  is  initially  zero.  Thus  the  total  number  of  such  pushes  is  at  most 
4A^i2  +  2n2  [log*  1/  +  lj  .□ 

As  in  the  Goldberg-Taijan  and  Ahuja-Orlin  algorithms,  the  time  to  perform  saturating 
pushes  and  relabeling  operations  is  O(nm).  The  only  significant  remaining  issue  is  how  to  choose 
vertices  to  add  to  S  when  S  is  empty.  For  this  purpose,  we  maintain  a  data  structure  consisting  of 
a  collection  of  doubly  linked  lists  list(r)  =  [i  e  N:  e(i)  >  A/k  and  d(i)  =  r]  for  each 
r  e  { l,2,...,2n-l }.  We  also  maintain  a  pointer  to  indicate  the  largest  index  r  for  which  list  ( r )  is 
non-empty.  Maintaining  this  structure  requires  0(1)  time  per  push  operation  plus  time  to  main¬ 
tain  the  pointer.  The  pointer  can  only  increase  due  to  a  relabeling  (by  at  most  the  amount  of 
change  in  label)  or  due  to  a  phase  change  (by  at  most  2 n).  Consequently,  the  number  of  times  the 
pointer  needs  to  be  incremented  or  decremented  is  0(n 2  +  nflog*  (/  +  !)).  The  overall  running 
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time  of  the  algorithm  is  thus  0(nm  +  kn2  +  n2(log*  U  +  1)).  Choosing  k  =f  log  U /log  log  t/] 
gives  the  following  result: 


Theorem  3.3.  The  stack-based  scaling  preflow  algorithm,  with  an  appropriate  choice  of  k,  runs  in 


0(nm  +  n 2 


log U 
log  logi/ 


)  time. 


Remark.  This  bound  can  be  improved  to  O  (run  +rt2  log £/*/loglogl/*),  where 
U*  -  4+  £  u  (s,v)/n.  The  algorithm  must  be  changed  slightly.  We  use  a  scaling  factor  of 

v:(i.v)e  E 

k  =  2  until  A<U*  and  then  switch  to  A:  log£/*/log  logi/*]  for  the  remainder  of  the  algo¬ 
rithm.  When  the  switch  occurs,  we  increase  A  (if  necessary)  to  be  a  power  of  k.  The  analysis  of 
the  initial  phases  (those  in  which  A  >  V*)  is  the  same  as  that  in  the  remark  following  Lemma  3. 1 . 
These  phases  account  for  only  O  ( n  2)  nonsaturating  pushes.  The  analysis  of  the  remaining  phases 
is  as  above  with  U  replaced  by  V* .  □ 


4.  Use  of  Dynamic  Trees 


The  approach  taken  in  Section  3  was  to  reduce  the  total  number  of  pushes  by  choosing  ver¬ 
tices  for  push/relabel  steps  carefully.  An  orthogonal  approach  is  to  reduce  the  total  time  of  the 
pushes  without  necessarily  reducing  their  number.  This  can  be  done  by  using  the  dynamic  tree 
data  structure  of  Sleator  and  Tarjan  [16,17,18],  We  conjecture  that,  given  a  version  of  the 
preflow  algorithm  with  a  bound  of  p>  mn  on  the  total  number  of  pushes,  the  running  time  can  be 


reduced  from  Ofp)  to  0(nm  log(-^-  +  1))  by  using  dynamic  trees.  Although  we  do  not  know 

mn 

how  to  prove  a  general  theorem  to  this  effect,  we  have  been  able  to  obtain  such  a  result  for  each 
version  of  the  preflow  algorithm  we  have  considered.  As  an  example,  the  0(nm  log  (n2/m)) 
time  bound  of  Goldberg  and  Tatjan  results  from  using  dynamic  trees  with  the  FIFO  selection 
rule;  the  bound  on  the  number  of  pushes  in  this  case  is  0(n3).  In  this  section  we  shall  show  that 
the  same  idea  applies  to  the  scaling  preflow  algorithm  of  Section  3,  resulting  in  a  time  bound  of 


nlogU 


+  2)). 


0(nm  log(- 

m  log  logl/ 

The  dynamic  tree  data  structure  allows  the  maintenance  of  a  collection  of  vertex-disjoint 
rooted  trees,  each  edge  of  which  has  an  associated  real  value.  Each  tree  edge  is  regarded  as  being 
directed  from  child  to  parent.  We  regard  every  vertex  as  being  both  an  ancestor  and  a  descendant 
of  itself.  The  data  structure  supports  the  following  seven  operations: 


find  root  (v):  Find  and  return  the  root  of  the  tree  containing  vertex  v. 


*-  <•-  a  s. 


find  size  (v): 


Find  and  return  the  number  of  vertices  in  the  tree  containing  vertex  v. 


find  value  (v):  Find  and  return  the  value  of  the  tree  edge  leaving  v.  If  v  is  a  tree  root, 

the  value  returned  is  infinity. 

find  min  (v):  Find  and  return  the  ancestor  w  of  v  with  minimum  find  value  (w).  In 

case  of  a  tie,  choose  the  vertex  w  closest  to  the  tree  root. 

change  value  (y,x):  Add  real  number  x  to  the  value  of  every  edge  along  the  path  from  v  to 

find  root  (v). 

link  (v,w,x):  Combine  the  trees  containing  v  and  w  by  making  w  the  parent  of  v  and 

giving  the  edge  (v.iv)  the  value  x.  This  operation  does  nothing  if  v  and 
w  are  in  the  same  tree  or  if  v  is  not  a  tree  root. 

cut  (v):  Break  the  tree  containing  v  into  two  trees  by  deleting  the  edge  joining  v 

to  its  parent;  return  the  value  of  the  deleted  edge.  This  operation 
breaks  no  edge  and  returns  infinity  if  v  is  a  tree  root. 

A  sequence  of  /  tree  operations,  starting  with  an  initial  collection  of  singleton  trees,  takes 
0(1  log  (z  +  1»  time  if  z  is  the  maximum  tree  size  [16,17,18]. 

In  our  application,  the  dynamic  tree  edges  form  a  subset  of  the  current  edges  out  of  all  the 
vertices.  Every  dynamic  tree  edge  is  an  eligible  edge;  its  value  is  its  residual  capacity.  The 
dynamic  tree  data  structure  allows  flow  to  be  moved  along  an  entire  path  at  once,  rather  than 
along  a  single  edge  at  a  time.  The  following  version  of  the  preflow  algorithm,  called  the  tree 
preflow  algorithm,  uses  this  ability.  Two  parameters  govern  the  behavior  of  the  algorithm,  a 
bound  A  on  the  maximum  excess  at  an  active  vertex,  and  a  bound  z,  1  <  z  <,  n,  on  the  maximum 
size  of  a  dynamic  tree.  The  algorithm  consists  of  beginning  with  an  initial  preflow  /  and  a  valid 
labeling  d,  and  with  every  vertex  as  a  one-vertex  dynamic  tree,  and  repeating  the  following  step 
until  there  are  no  active  vertices: 

tree  pushJrelabel  (v): 


Applicability:  Vertex  v  is  active. 


Action:  Let  (v,w)  be  the  current  edge  out  of  v.  Apply  the  appropriate  one  of  the 
following  cases: 

Case  1:  (v,w)  is  eligible.  Let  x  =  find  root  (w),  e  =  min  (e(v),  uf  ( v,w),find  min  (w)), 
and  8  =  e  if  *  =  r,  8  =  min  {e,  A  -e(x)}  if  x  *  t.  Send  8  units  of  flow  from  v  to  x  by 
increasing  f(v,w )  by  8  and  performing  change  value  (w,-8).  This  is  called  a  tree  push 
from  v  to  x.  The  tree  push  is  saturating  if  8  =  min  {  Uf  (v,w),find  min  (v)}  (before  the 
push)  and  nonsaturating  otherwise.  If  the  tree  push  is  nonsaturating,  e(v)  =  0  (after  the 
push),  and  find  size  (v)  +  find  size  (w)  £  z,  perform  link  ( v,w,Uf  (v,w)).  Otherwise,  while 
find  value  ( find  min  (w))  =  0  perform  cut  (find  min  (w)). 

Case  2:  (v,w)  is  not  eligible  and  not  last  on  A(v).  Replace  (v,w)  as  the  current  edge 
out  of  v  by  the  next  edge  on  A(v). 

Case  3:  (v.w)  is  not  eligible  and  is  last  on  A(v).  Relabel  v.  Replace  (v,w)  as  the 
current  edge  out  of  v  by  the  first  edge  on  A(v).  For  every  tree  edge  (y,  v)  perform 
cut  (y). 

This  algorithm  stores  flow  in  two  different  ways:  explicitly  for  edges  that  are  not  dynamic 
tree  edges  and  implicitly  in  the  dynamic  tree  data  structure  for  edges  that  are  dynamic  tree  edges. 
After  each  cut,  the  flow  on  the  edge  cut  must  be  restored  to  its  correct  current  value.  In  addition, 
when  the  algorithm  terminates,  the  correct  flow  value  on  each  remaining  tree  edge  must  be  com¬ 
puted.  For  edges  cut  during  the  computation,  the  desired  flow  values  are  easily  computed,  since 
the  current  residual  capacity  of  an  edge  that  is  cut  is  returned  by  the  cut  operation.  Computing 
correct  flows  on  termination  can  be  done  using  at  most  n  find  value  operations. 

The  algorithm  just  presented  is  a  generic,  modified  version  of  the  dynamic  tree  algorithm  of 
Goldberg  and  Taijan  [8,9].  It  maintains  the  following  invariants:  every  active  vertex  is  a  tree 
root;  no  excess  exceeds  A;  no  tree  size  exceeds  z;  every  tree  edge  is  eligible.  The  analysis  of 
Goldberg  and  Taijan  gives  the  following  result: 

Lemma  4.1  [9].  The  total  time  required  by  the  dynamic  tree  preflow  algorithm  is 
0(nm  log  (z  +  1))  plus  0(log  (z  +  1))  per  tree  push.  The  number  of  links,  cuts,  and  saturating 
tree  pushes  is  O(nm). 

The  efficiency  of  the  algorithm  depends  on  the  number  of  nonsaturating  tree  pushes,  which 
in  turn  depends  on  how  tree  push/relabel  steps  are  selected.  We  shall  show  that  a  selection  order 

n2 

analogous  to  that  used  in  Section  3  results  in  a  bound  of  0(nm  +kn2  +  —  (log*  U  +  1))  on  the 


* 


number  of  nonsaturating  tree  pushes,  where  A:  is  the  scaling  factor.  From  this  a  total  time  bound 

of  0(nm  log  ( — 7  +  2))  follows  by  an  appropriate  choice  of  k  and  z. 

m  log  logt/ 

We  call  the  resulting  version  of  the  preflow  algorithm  the  treelscaling  preflow  algorithm. 
In  addition  to  parameters  A  (the  maximum  excess)  and  z  (the  maximum  tree  size),  the  algorithm 
uses  a  third  parameter  k,  the  scaling  factor,  and  a  stack.  S.  Whereas  z  and  k  remain  fixed 
throughout  the  running  of  the  algorithm,  A  decreases.  Initially  S  =  $  and  A  is  the  smallest  integer 
power  of  k  such  that  A  £  U.  An  active  vertex  v  is  said  to  have  large  excess  if  e(y)  >  A  Ik  and 
small  excess  otherwise.  A  dynamic  tree  is  said  to  be  large  if  it  contains  more  than  z/2  vertices 
and  small  otherwise.  The  tree/scaling  algorithm  consists  of  repeatedly  applying  the  tree 
push/relabel  step  to  the  top  vertex  on  S,  while  manipulating  S  according  to  the  following  rules: 

Rule  1.  If  5  is  empty,  add  to  S  any  active  vertex  that  is  the  root  of  a  small  tree.  If  all  active  ver¬ 
tices  are  roots  of  large  trees,  add  to  S  a  large-excess  active  vertex  of  largest  label.  If  all  active 
vertices  are  roots  of  large  trees  and  all  have  small  excess,  replace  A  by  A  Ik.  (The  algorithm  ter¬ 
minates  when  A  <  1.) 

Rule  2.  After  a  tree  push  from  a  vertex  v  to  a  vertex  x,  if  the  push  is  nonsaturating  and  e(v)  >  0, 
push  x  onto  S.  Otherwise,  while  S  is  nonempty  and  its  top  vertex  has  small  excess,  pop  S. 

Rule  3.  After  relabeling  a  vertex,  pop  it  from  S,  and  continue  to  pop  S  until  it  is  nonempty  or  its 
top  vertex  has  large  excess. 

Every  tree  push  made  by  the  tree/scaling  algorithm  is  from  a  vertex  that  is  either  the  root  of 
a  small  tree  or  has  large  excess.  As  discussed  in  Section  3,  the  time  necessary  for  implementing 
Rules  1-3  is  0(n2  +  nlog*  (U  +  1))  plus  0(1)  per  tree  push,  which  is  dominated  by  the  bound  on 
tree  pushes  derived  below. 

2  n2 

Lemma  4.2.  The  tree/scaling  algorithm  makes  0(nm  +  kn  +  —  Gog*  U  +  1))  nonsaturating 
tree  pushes. 

Proof.  The  argument  is  analogous  to  the  proof  of  Lemma  3.2.  Let  the  potential  O  be 
defined  by  =  £  «00  <*(v)/A.  We  have  0  £  <  n2.  Increases  in  <I>  are  due  to  relabelings 

v  active 

and  changes  in  the  value  of  A.  The  total  increase  in  $  due  to  relabelings  is  0(n2).  Immediately 
after  a  change  in  the  value  of  A,  O  S  2nz/z,  because  every  active  vertex  must  be  the  root  of  a 
large  tree,  and  only  2/i/z  large  trees  can  coexist.  Thus  the  total  increase  in  4>  caused  by  changes 


in  A  is  0(—  Gog*  £/  + 1)). 

Consider  nonsaturating  tree  pushes  that  add  a  vertex  to  S.  For  such  a  push,  let  x  be  the  ver¬ 
tex  added  to  S.  We  charge  the  push  either  to  the  relabeling  of  x  that  pops  x  from  S  or  to  the  at 
least  1-1  Ik  decrease  in  <t>  caused  by  pushes  from  x  before  x  is  popped.  The  total  number  of  such 
n2 

pushes  is  thus  0(n2  +  —  (log*  U  +  1)).  Similarly,  the  number  of  nonsaturating  tree  pushes  that 

begin  with  more  than  one  vertex  on  S  and  pop  at  least  one  vertex  from  S  is  also 
_2 

0(n2  +  Gog*  U  +  1)). 

The  last  kind  of  nonsaturating  tree  push  that  we  must  count  is  a  push  from  a  vertex  v  that  is 
the  only  vertex  on  S  and  is  popped  from  S  after  the  push.  We  call  such  a  push  an  emptying  push. 
An  emptying  push  from  a  vertex  v  reduces  e(v)  to  zero.  To  count  emptying  pushes,  we  divide 
them  into  those  from  roots  of  small  trees  and  those  from  roots  of  large  trees. 

We  charge  an  emptying  push  from  the  root  v  of  a  small  tree  to  the  event  that  either  created 
the  tree  or  last  made  v  active,  whichever  occurred  first.  If  the  tree  was  created  after  v  became 
active,  we  charge  the  push  to  the  link  or  cut  that  created  the  tree.  Each  link  or  cut  is  charged  to  at 
most  one  push,  and  there  are  O(nm)  links  and  cuts,  which  gives  an  0(nm )  bound  on  such  pushes. 
If  the  tree  was  created  before  v  became  active,  then  v  became  active  as  a  result  of  a  saturating 
push,  of  which  there  are  O(nm),  or  as  a  result  of  a  nonsaturating  push  from  the  root  of  a  large 
tree.  We  have  already  counted  the  number  of  such  pushes  that  are  nonemptying.  Thus  it  suffices 
to  count  emptying  pushes  from  the  roots  of  large  trees. 

We  count  such  pushes  using  a  second  potential  <X>2-  To  define  d>2.  we  need  to  define  an 
integer  /  and  a  set  P,  as  in  the  proof  of  Lemma  3.2.  The  value  of  /  is  the  minimum  of  2 n  and  the 
smallest  label  of  a  large-excess  root  of  a  large  tree  added  to  S  while  S  was  empty  during  the 
current  phase.  (As  in  Section  3,  we  define  a  phase  to  consist  of  a  maximal  period  of  time  during 
which  A  stays  constant.)  The  set  P  consists  of  all  vertices  with  label  greater  than  /  and  all  vertices 
with  label  equal  to  /  from  which  an  emptying  push  has  been  made  during  the  current  phase.  The 
definition  of  4%  is 

<X>2  =  £  e(v)  (d(v)  -/  +  1)/A.  (lfP=0,Q>2=  0  ) 

* *  p  ■■ «(»)  >  o 


The  total  increase  in  <r>2  over  all  phases  is  0(n2  +  —  Gog*  U  +  1)).  The  0(n2)  contribu- 

n 2 

tion  is  due  to  relabelings.  The  0(—  Gog*  U  + 1))  contribution  is  due  to  decreases  in  l.  A 


decrease  in  l  by  one  causes  an  <9(-^j)  increase  in  d>2,  since,  just  after  such  a  decrease  in  /  occurs, 

every  vertex  in  P  has  level  greater  than  l,  and  every  active  vertex  in  P  has  small  excess  and  is  the 
root  of  a  large  tree. 

Every  emptying  push  from  the  root  of  a  large  tree  causes  either  a  decrease  of  at  least  1  Ik  in 
<b2,  which  can  happen  0(/fc/r  +  —  (log*  U  +  1»  times,  or  an  addition  of  a  vertex  to  P,  which 
can  happen  0(n(\ogt  U  +  1))  times.  Since  z  <.  n,  the  latter  bound  is  dominated  by  the  former. 
Combining  all  our  estimates  gives  the  lemma.  □ 

kn2 

Remark.  A  slightly  weaker  bound  of  0(nm  +  kn2  + - (log2  t/  +  1))  on  the  number  of  nonsa¬ 

turating  tree  pushes  can  be  obtained  with  a  simpler  argument  using  only  <I>,  analogous  to  the 
proof  of  Lemma  3.1.  □ 


Theorem  4.3.  With  a  choice  of  k  ~\  —  +  2]  and  z  =  min  {["  —  log*  f/  + 1  j  ,n)  the  tree/scaling 

n  m 

preflow  algorithm  runs  in  0(nm\og  ( — ”  *°j^.  +  2))  time. 

mloglogU 

Proof.  By  Lemmas  4.1  and  4.2,  the  running  time  of  the  algorithm  is 
2 

0{(nm  +  kn2  +  —  Gog*  U  +  1))  log(z  +  1)).  This  is  0(nm  log n)  if  z  =  [  —  log*  U  +  l"|  <  n. 
i  m 

On  the  other  hand,  if  z  =  n  the  computation  is  complete  after  one  phase,  and  the  running  time  of 

the  algorithm  is  again  0(nm  logn).  We  consider  two  cases. 

Case  1:  logf/ /  log  logi/ >  m/Vn.  Then  log( — ^ -  +  2)  =  Cl  Oogn),  and  the 

mloglogC/ 

theorem  is  true. 

Case  2:  logf/  /  log  logf/  S  ml'fn.  Then  logjtZ/k^logjf/  S«3/<,  which  implies 

nl 

log2l/  =  0(m).  Let  us  estimate  the  quantity  —  Gog*  V  + 1).  If  z  =  n,  this  quantity  is 

n  Gog*  U  +  1)  =  0(nm).  On  the  other  hand,  if  z  =  [  —  log*  U  +  l"|  , 

m 

n2 

—  Gog*  U  + 1)  £  2nm.  It  follows  that  the  running  time  of  the  algorithm  is 

nlog*(/ 

0(nmlog  (z  +  l))  =  0(nmlog( - +  2)).  It  remains  to  prove  that 

m 

log  (  Wl0g—  +  2)  =  OQog  (—7^77  +  2)).  We  have  log*  U  =  log U  /  log  (—  +  2). 
m  m  loglog U  n 


Again  the  proof  divides  into  two  cases: 


Case  2.1:  log U  >  (—  +  2)4.  Then  -j-  log  logf/  >  log  (—  +  2)  >  log (m/n).  Also, 
b  4  n 

log  logU  >  4,  which  implies  log  log  log U  <  y  log  logt/.  Hence 

i°g  (  7 +  2)  S  log  log U  -  log  log  log U  -  log  (m/n)  £  —  log  log U.  But 

mlog  logt/  4 

then  log  (— +  2)  <  log  (logt/  +  2)  =  <9(log( — "  — ^  +  2)),  as  desired. 
m  cc  m  log  logt/ 


Case  2.2:  logt/  £  ( —  +  2)4.  In  this  case  log  logt/  <,  4  log  (—  +  2),  which  implies 

log  + 2) = log  + 2)  £  log  (.  +  2) , 

m  m  .  m  log  log!/  m  log  logt/ 

m  log(—  + 1) 

as  desired. 


Thus  in  all  cases  the  algorithm  runs  in  the  claimed  time  bound.  □ 

Remark.  The  bound  in  the  theorem  can  be  improved  to  (9  (run  log  ( — n  — -  +  2),  analo- 

m  log  logt/ 

gously  to  the  improvements  in  the  bounds  of  Section  3.  To  obtain  this  improvement,  we  run  one 
of  the  algorithms  of  Section  3  with  k  =  2  until  A  <  t/*,  and  then  switch  to  the  dynamic  tree  algo¬ 
rithm.  □ 

5.  Remarks 

The  scaling  preflow  algorithm  of  Section  3  not  only  has  a  good  theoretical  time  bound,  but 
is  likely  to  be  very  efficient  in  practice.  We  intend  to  conduct  experiments  to  determine  whether 
this  is  true.  The  algorithm  of  Section  4  is  perhaps  mainly  of  theoretical  interest,  although  for 
huge  networks  there  is  some  chance  that  the  use  of  dynamic  trees  may  be  practical. 

The  two  obvious  open  theoretical  questions  are  whether  further  improvement  in  the  time 
bound  for  the  maximum  flow  problem  is  possible  and  whether  our  ideas  extend  to  the  minimum 
cost  flow  problem.  It  is  not  unreasonable  to  hope  for  a  bound  of  0(nm )  for  the  maximum  flow 
problem;  note  that  the  bounds  of  Theorems  3.3  and  4.3  are  0(nm )  for  graphs  that  are  not  too 

sparse  and  whose  edge  capacities  are  not  too  large,  i.e.  =  0(mln ).  Obtaining  a  bound 

better  than  0(nm)  would  seem  to  require  major  new  ideas. 


As  a  partial  answer  to  the  second  question,  we  have  been  able  to  obtain  a  time  bound  of 
0(nm  loglog  U  log  (nC))  for  the  minimum  cost  flow  problem,  where  C  is  the  maximum  edge 
cost,  assuming  all  edge  costs  arc  integral.  This  result  uses  ideas  in  the  present  paper  and  some 
additional  ones.  It  will  appear  in  a  future  paper. 
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